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ABSTRACT 

Kepler's supernova, discovered in October 1604, produced a remnant that has been well studied 
observationally in the radio, infrared, optical, and X-ray bands, and theoretically. Some models have 
predicted a TeV gamma-ray flux that is detectable with current Imaging Cherenkov Atmospheric 
Telescopes. We report on observations carried out in 2005 April with the CANGAROO-III telescope. 
No statistically significant excess was observed, and limitations on the allowed parameter range in the 
model are discussed. 

Subject headings: gamma rays: observation — supernova: individual (Kepler's SNR) 



1. INTRODUCTION 

Kepler's supernova remnant (SNR) (G4.5+6.8) is 400 
years old (pain [20051 see for review) and provides an 
unrivaled opportunity to verify the belief that super- 
nova remnants are the origin of Galactic cosmic rays. 
Cas A is younger (b y ~60 years) and was detected at 
TeV 7-ray energies (EE aroman et al.ll2001h . which im- 
plies the acceleration of high-energ y cosmic rays. Older 
remnants, such as RX J0852.0-4622 (iKatagiri et al.ll2005l : 
lAharonian et alj 120051: lEnomoto et all l2006bl) and RX 



J171 3. 7-3946 (lEnomoto et all l2002aMAharonian et all 
120041 ) . both thought to be 1,000~2,000 years old, have 
also been detected. If SNR age was the dominant factor 
in for cosmic-ray acceleration, one might expect simi- 
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lar levels of cosmic ray acceleration in Kepler's SNR. Of 
course, other variables, such as SN type, local environ- 
ment, and distance, will certainly have some impact on 
the likelihood of TeV gamma-ray detection from a SNR. 

Kepler's SN was considered to be a type la super- 
nova (SN) based on an inter pretation of the histor- 
ical light curve (iBaadd [1943). It was later shown 
that the light cu rve was also in agreem ent with 
a type II-L SN (jDoggett k Branchil I1985D . How- 
ever, recen t observations of thermal X-ray emission 
by ASCA (iKimigasa k TsunemU I2007D and Chandra 
([Reynolds et all 120071) demonstrated that the SNR re- 
sulted from a thermonuclear supernova (type la), rather 
than the core-collapse of massive star (type II), even 
though there is evidence that the remnant is interacting 
with the progenitor star wind material. It may be that 
a type la event took pla ce in a more massive progenitor 
st ar with a strong wind (|Revnolds et al.ll2007D . 

iBerezhko et al.1 (|2006h have modeled Kepler's SNR and 
predicted a detectable TeV gamma-ray flux under vari- 
ous assumptions on distances and supernova kinetic ener- 
gies that had been previously discussed in the literature. 
Their prediction can be probed by the H.E.S.S. 17 and 
the future GLAST 18 experiments. The CANGAROO- 
III imaging Cherenkov atmospheric telescope is less sen- 
sitive by a factor of 3-5 than H.E.S.S.; however, it is 
able to study a specific parameter range of the models, 
i.e., around the region of a supernova explosion energy 
of 10 51 erg and distance of 4.8 kpc. Here, we report on 
the re sult of the 2005 April observations. As extensions 
to the Bcrczhk o et al.l (|2006l ) theory have successfully ex- 
plained the gamma-ray fluxes from other historical SNRs, 
it is important to investigate and constrain the allowed 
parameter ranges in this model with measurements of 
fluxes or upper limits. 



2. CANGAROO-III STEREOSCOPIC SYSTEM 
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18 See 
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The CANGAROO-III stereoscopic system consists of 
four imaging atmospheric Cherenkov telescopes located 
near Woomera, South Australia (31°S, 137°E). Each 
telescope has a 10 m diameter segmented reflector, con- 
sisting of 114 spherical mirro rs made of fiber-reinforced 
plastic (jKawachi et alJ l2001h , each of 80 cm diameter, 
mounted on a parabolic frame with a focal length of 
8 m. The total light-collecting area is 57.3 m 2 . The first 
telescope, Tl, which w as the CANGAROO-II telescope 
(jEnomoto et al.|[2002a( ) , is not presently in use due to its 
smaller field of view and higher energy threshold. The 
second, third, and fourth telescopes (T2, T3, and T4) 
were operated for the observations described here. The 
camera syst ems for T2, T3 , and T4 are identical and are 
described in lKabuki et al. (2003). The telescopes are lo- 
cated at the eastern (Tl), western (T2), southern (T3) 
and nor thern (T4) corners of a diamond with sides of 
~100m (jEnomoto et al.ll2002bf ). The point-spread func- 
tions of these telescopes are 0.°24. 

3. OBSERVATIONS 

The observations were carried out during the period 
from 2005 April 11 to 17 (UT) using the "wobble mode" 
in which the pointing position of each telescope was 
shifted in declinat ion by ±0.5 degree every 20 minutes 
(jDaum et al.lll997ft from the target: (RA, dec [J2000]) = 
(262.°671, -21.°486). We made no OFF source runs, as 
the wobble mode enables OFF-source regions to be ob- 
served simultaneously with the target regions. The sen- 
sitive region in wobble mode observations is considered 
as being within one degree from the average pointing po- 
sition. This SNR is located 6.°8 from the Galactic plane; 
therefore, no significant diffuse gamma-ray background 
is expected within the field of view. 

In the observations, the hardware trig ger used to select 
any two telescope hits was employed ([Nishiiima et al] 
2005). The images in two out of three telescopes were 
required to have clusters of at least five adjacent pixels 
exceeding a 5 photoelectron threshold (off-line two-fold 
coincidence). To illustrate the effect of this criterion, the 
event rate was reduced from 10~12 to 6^8 Hz for T3-T4 
coincidences, depending on the elevation angle. Looking 
at the time dependence of these rates, we can remove 
data taken under cloudy conditions. The effective obser- 
vation time was 874 minutes, and the mean zenith angle 
was 15. °2. 

The light-collecting efficiencies, including the reflec- 
tivity of the segmented mirrors, the light guides, and 
the quantum efficiencies of the photo multiplier tubes 
were m onitored by a muon-ring analysis (jEnomoto et al.l 
12006a!) with individual trigger data during the same pe- 
riod. The average light yield per unit arc-length of muon 
rings is approximately proportional to the light-collecting 
efficiencies. Deterioration in these efficiencies is mostly 
due to dirt and dust settling on the mirrors and light 
guides, which are washed annually to improve their re- 
flectivities. In analyzing T2 data, we had some difficul- 
ties in detecting muon-rings during this period; there- 
fore, we did not use T2 in this analysis. Unfortunately 
these observations were made shortly before regular mir- 
ror washing. Also we had some mis-setting of the T2 
ADC-gate width in this period. This analysis, therefore, 
used only T3 and T4 two-fold coincidence data. 



4. ANALYSIS 

The analysis pr ocedures used here w ere identical to 
those described in Ka buki et al.l (|2007[ ). except for the 
point that these were two -fold coincidenc e data. More 
details can be found in I Enomoto et alJ (|2006aD and 
lEnomoto et al . ( 2006b). Here, we briefly describe them. 

At first, t he image moments of Width and Length 
(|Hillaslll985n were calculated for the two telescopes. The 
incident direction of the gamma-ray was determined by 
minimizing the sum of the squared widths (weighted by 
the photon yield) of the two images seen from the as- 
sumed position (fitting parameter) with a constraint on 
the distances from the intersection point to each image 
center. 

In order to derive the gamma-ray likeli ness, we used 
the Fisher Discrimina nt (hereafter FD) ((Fisher 1936; 
lEnomoto et al.l[20 06a). The input parameters were 

P= (W3,W4,L3,L4), 

where W3, WA, L3, LA are energy-corrected Widths and 
Lengths for T3 and T4. 

We rejected events with any hits in the outermost layer 
of the cameras ("edge cut"). These rejected events were 
potentially incompletely sampled, resulting in errors par- 
ticularly in the Length distribution, which would have 
produced deformations of the FD. 

Then FD distributions were derived on a position-by- 
position basis. Comparing those in the signal region and 
the control background region with the Monte-Carlo ex- 
pectation, we can derive the number of gamma-ray-like 
events. Here, we assume the FD distribution of the 
gamma-ray signal to be that derived from Monte-Carlo 
simulations. In the gamma-ray simulations we used a 
spectrum proportional to E 1 , where 7=— 2.1±0.2. Fits 
of the FD distribution of the source position with the 
above simulated signal and control background functions 
were carried out to derive the number of gamma-ray-like 
events. This was a one-parameter fitting with the con- 
straint that the sum of the signal and the background 
events corresponds to the total number of events, i.e., 
the fitted parameter can be derived exactly analytically. 

5. RESULTS 

Since the spatial size of Kepler's SNR (100") is much 
less than our angular resolution (a = 0.24°), we concen- 
trate here on searching for a point source near the target 
center. 

In order to determine whether or not there is a gamma- 
ray excess around the SNR, we made the FD distribu- 
tion within the PSF (6> 2 < 0.06 ~ 0.24 2 ) and fitted it 
with a background function derived from the FD distri- 
bution in the region 9 2 = (0.1-0.2), and a signal function 
from Monte-Carlo simulations. The fitting parameter 
is the ratio of gamma-rays to total number of events. 
The fitting results are shown in Fig. [T] The best-fit 
excess was 71±32 events, where the uncertainty is the 
la statistical error. The gamma-ray response function 
from the Monte-Carlo simulation is shown by the red 
histogram. The threshold of this analysis is estimated 
from the Monte-Carlo simulation to be ~500GeV. The 
systematic error on the energy determination is consid- 
ered to be less than 15%, with the main factor being the 
uncertainty in the light collection efficiency and atmo- 
spheric conditions. 
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Fig. 1. — FD distribution for the region inside 2 < 0.06 dog 2 . 
The black points with error bars are those for the above region, 
the green histogram is for e 2 inside the (0.1-0.2) [deg 2 ] region, the 
blue points with error bars are subtracted data using the results of 
the fit, that is, "gamma-ray-like" events, and the red histogram is 
the best-fit gamma-ray, i.e., the gamma ray response function. 



S; SO 



60 -; 



40 :};-■ 



20 



-20 - 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

6" (degree ) 

Fig. 2. — 9 2 plot in the unit of deg 2 . The points with error 
bars are the result of fit, that is, the "gamma-ray-like" event dis- 
tribution as a function of 9 2 . The dot-dashed line is the zero level. 
The (light) dotted histogram is the best-fit for our point-spread 
function. The (heavy) dashed histogram is the two-er upper limit 
for the point-source assumption. 

We then made a radial distribution of gamma-ray-like 
events. FD distributions in various 9 2 slices were made. 
The control background region was again selected in the 
8 2 range of between 0.1 and 0.2 deg 2 . The standard 
fitting procedure described above was carried out. The 
fitted result is shown in Fig. [2] The reduced \ 2 f° r a nuu 
assumption (the dot-dashed line) is x 2 /DOF = 25.4/25 
(where DOF is degrees of freedom). The best fit with 
the point spread function (PSF) is shown by the dashed 
histogram, where x 2 /D.O.F = 18.1/25. The dashed his- 
togram is the 2cr upper limit (135-event excess) for the 
PSF excess. 

In order to examine the morphology, we segmented 
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Fig. 3. — Significance map. The average telescope pointing 
position is indicated by the white cross at the center. The dotted- 
white circle is the point-spread function and the searched region. 
The dashed- white circle is the fiducial region (1° radius). The thin- 
white contours are the radio measurement at 4850 MHz ()SkvViewl 
l2007t ) which are well inside the searched region. The color bar 
indicates the excess in standard deviations over the background 
(see text for details). 

the field of view into 0.2 x 0.2 degree 2 square bins. The 
FD distributions for corresponding bins were made and 
fitted. The control-background region is defined as the 
second-closest layer of 16 bins, all of which are more than 
0.3 deg from the center of the target region, i.e., larger 
than the 0.24 degree point-spread function (PSF). The 
statistics of the control-background are, therefore, six- 
teen times larger than that of the signal bin. The results 
are shown in Fig. [3j We smoothed the results by aver- 
aging the neighboring nine bins. Our sensitivity falls off 
significantly beyond one degree in radius from the cen- 
ter. The PSF is the 0.24 degree radius circle which fully 
contains the radio observed SNR (the thin-white con- 
tours). According to the Monte-Carlo simulations, 65% 
of gamma-rays from this SNR should be contained in 
this circle. The PSF is not a Gaussian and has a broader 
tail component. In order to contain 90% of events, we 
need to broaden this cut to 0.5 degree, resulting in a 
loss of sensitivity. We, therefore, selected the cut at 0.24 
degree (la region). The significance distributions (ex- 
cess divided by the statistical error before smoothing) 
are approximately normal (Gaussian) distributions with 
a mean value of 0.21±0.13, and a standard deviation of 
1.19±0.12, consistent with null assumption within sys- 
tematic uncertainties. The statistical significance of the 
maximum located 0.°35 west-north-west from the center 
is (before smoothing) 3.3<t, and therefore not compelling. 

There is no statistically-significant excess that is con- 
sistent with emission from Kepler's SNR convolved with 
the telescope PSF. The derived upper limits (ULs) for 
the gamma-ray flux are shown in Table [1] Here, we used 
a E~ 2A±a - 2 spectrum for the gamma-ray simulation. The 
ULs range between 10-30% of the Crab nebula flux. The 
statistically insignificant excess near the center of the 
field of view, shown in Figs.HJEJ and [3] only appeared in 
the lower energy regions. At higher energies, we do not 
see any excess. Therefore the ULs at lower energies were 
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TABLE 1 

The 2a upper limits to the integral fluxes from Kepler's 
snr at five energy thresholds. the spectral index, 7, of 
the energy spectrum (e _7 ) is that used in the 
Monte-Carlo gamma-ray simulations. Note that the 
gamma-ray acceptance depends on 7. 



7 Threshold [GcV] Upper Limits [cm 2 s *j 
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Fig. 4. — Spectral energy distributions for a fixed distance of 
4.8 kpc. The black curves with arrows at both ends are the upper 
limits obtained by these observations. Three curves are obtained by 
different acceptance corrections using three types of energy spectra 
in the Monte-Carlo simulation (see the caption of Table [TJ. The 
colored cur ves are the theoretical predictions, which are as same 
as Fig. 3 of lBerezhko et al.1 J20TJ6T) . 

higher than that at higher energies. 

6. DISCUSSION 

The upper limits given in Table Q] are plotted on spec- 
tral energy distributions in Figs. H] and O The verti- 
cal and horizontal un its were fitted to Figs. 3 and 4 in 
iBerezhko et al.l (|2006f) in order to discuss the allowed pa- 
rameter ranges with respect to our observational upper 
limits. This theory considered a reasonably wide range 
of possibilities for the distance of this SNR (3.4-7 kpc) 
and the supernova explosion energy (0.5-2 x 10 51 erg). 
Other adopted parameters, such as the cosmic ray in- 
jection rate, expansion rate, and electron-to-proton ra- 
tio, while plausible, are open to debate. We do not 
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Fig. 5. — Spectral energy distributions for a fixed supernova 
explosion energy of 10 51 erg. The black curve with arrows at both 
ends are the upper limits obtained by this observation. Three 
curves are obtained by different acceptance corrections using three 
types of energy spectra in the Monte-Carlo simulation (Table [TJ. 
The colored cur ves are the theo r etical predictions, which are as 
same as Fig. 4 of Bcrczhko et al. (200(1). 

review th ose details here but r efer readers to the dis- 
cussion m IBerezhko et~all (|2006[ ). The black curve with 
arrows at both ends were obtained from this observa- 
tion. Three types of the energy spectra (oc E 1 : 7 = 
— 1.9, — 2.1, — 2.3) were used for acceptance correc- 
tion using the Monte-Carlo simulation. The uncertainty 
due to the assumption of the energy spectral index is 
small on a logarithmic scale. The colo red curves were 
obtain ed from a theoretical prediction by Bcrczhko et al.l 
(2006). In Fig. [4] the distance to the object was fixed to 
be 4.8 kpc. The red, green, blue, and yellow curves cor- 
respond to different supernova explosion kinetic energies 
and corresponding number densities of ambient circum- 
stellar material, which came from fitting to the observed 
shock radius and speed. The most probable is the sec- 
ond one (the green curve) and our upper limits are (in 
part) below it, meaning that for a distance of 4.8 kpc 
the explosion energy should be less than ~ 10 51 erg. Al- 
though the green curve is close to a best estimation, a 
large allowable range of parameter space remains. 

Fig. [5] is the case when the supernova explosion energy 
is fixed at 10 51 erg, and severa l distances in the allowed 
range dRevnoso fc Gossl [l999) are assumed. Distances 
less than ~5kpc are (for this fixed SN energy) not fa- 
vored, suggesting that this SNR is marginally more dis- 
tant than the best current observational estimations. Of 
course, this conclusion also assumes that all other as- 
sumptions in the model are correct, for example, that 
the expansion is in Sedov phase, that the supernova was 
Type la, that 10% of kinetic energy is transferred to the 
cosmic-ray energy, etc. 

The constraints on the parameters of the distance d 
and the ambient density npj are illustrated in Fig. [6l 
Here we assume that the total number spectrum of pro- 
tons is proportional to a power law with an exponential 
cutoff E~ p exp(— E/E max ) and that the neutral pion de- 
cay 7-ray emission dominates. If the conversion efficiency 
(e) from the explosion energy E sn to the cosmic-ray en- 
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Fig. 6. — The allowed region in the plane of the distance (d) 
versus ambient density (nn), based on the neutral pion decay 7-ray 
emission model where we assume that the total number spectrum 
of protons is proportional to a power law with an exponential cut- 
off and the value of the conversion efficiency from the explosion 
energy to the cosmic-ray energy is 10%. The solid, dashed, and 
dotted lines indicate the upper limits obtained from the observed 
2<r upper limits in TablelTI compared with the integral fluxes of the 
models for the power-law indices 2.1, 1.9, and 2.3, respectively. The 
colors of the lines represent the cutoff energy of protons; the black, 
red, and blue ones are obtained with the values 10 15 , 10 14 , and 
10 13 eV, respectively. The allowed regions lie under these upper- 
limit lines. The hatched region indicates the region which satisfies 
the Sedov- Taylor solution of the apparent radius of the Kepler's 
SNR, 100" at 400 yr with the explosion energy 0.5 ~ 2 X 10 s erg. 
The dot-dashed lines are obtained from the approximate analytic 
blast-wave positions of th e radius 100", assuming the solution of 
ITruelove fc McKeei l Tl999T) . The red and blue ones are obtained as- 
suming the ejecta power-law index n = 6 and n = 14, respectively, 
with ejecta mass of 1.4Mq with the explosion energy 10 erg. 

ergy is assumed to be 10% (i.e., e = 0.1), the normal- 
ization factor of the proton spectrum can be determ ined 
and the 7-ray fluxes can be calculated (|Moril [l997h on 
the assumption of the power-law index p and the cutoff 
energy E max . Given the parameters of the E sni p, and 
E max , the upper limits oinyi/d 2 are calculated from the 
observed 2a upper limits in Table [TJ compared with the 
integral fluxes of the model, because the 7-ray fluxes F 1 
are proportional to eE sn nn/d 2 . In Fig. we plotted 
the upper limits for the power-law indices 2.1 (solid), 1.9 
(dashed), and 2.3 (dotted) and for the cutoff energies of 
protons 10 15 (black), 10 14 (red), and 10 13 eV (blue). We 
did not plot the case of p — 2.3 and E max = 10 13 eV, as 
the flux in the GeV energy region exceeds the EGRET 
upper limit. 

The apparent radius 0=100" of the Kepler's SNR at 
i ag e gives another constraint. Here we consider two so- 
lutions on the expansion law of the blast-wave shock. 
The first one is the Sedov- Taylor solution : 9d oc 

(£ SI1 /n H ) 1/5 ^e . Another one is the ap proximate ana- 
lytic solutions (jTruelove &; McKeei fT999). which can be 
applied from the ejecta-dominated phase to the Sedov- 
Taylor phase. In the latter case, the extra parameters of 
ejecta mass M e j and the ejecta power-law index n, are 
added to the explosion energy E sn and the ambient mat- 
ter number density nn . The region which is satisfied with 
the Sedov- Taylor solution of the apparent radius of the 



Kepler's SNR, 100" at 400 yr and E sn = 0.5 - 2 x 10 51 
erg are shown as the hatched one in Fig. 6. The dot- 
dashed lines are obtained from the approximate analytic 
solution with ejecta mass of 1.4M Q with the explosion 
energy 10 51 erg with two kinds of the ejecta power-law 
index n — 6 and n = 14. In the case of p = 2.1 and 
E max = 10 14 eV, we note that the CANGAROO-III up- 
per limit implies that the Kepler's SNR is located at a 
distance larger than about 4.5 kpc. For the maximum 
upper limit in the case of p = 2.1 and E max — 10 13 eV, 
this means that the SNR is located at a distance larger 
than about 3.6 kpc. 

The distance of a type la supernova can be estimated 
using the correlation between the shape of the optical 
light curves and the intrinsic luminosity of SNe. Before 
knowledge of this correlation, iBaadd (|1943| ) studied the 
historical light curv e of the Kepler's S N an d classified 
it as a type la, and iDanziger &: Gossl (|1980l ) estimated 
the distance of 3.2 ± 0.7 kpc using only the maximum 
luminosity. We can now fit Baade's d ata with the im - 
proved light-curve template of a type la (|Jha et al.ll2007h : 
d = 4.0±0.4 kpc on the assump tion of a visual extinction 
3.27±0.14 mag (|Schaeferlll996h . although the fitted light 
curve after 100 days is not a good fit. 

This value seems to be marginally consistent with 
our lower limit obtained by this TeV 7-ray observa- 
tions. On the other hand, based on the study of Hi 
kine matics and the as s ociati on of Hi cloud with the 
SNR, iRevnoso fc Gossl (|1999h put a lower limit of d = 
4.8 ± 1.4 kpc and an upper limit of 6.4 kpc on the dis- 
tance. These estimations do not contradict our lower 
limit. 

In any case, a part of the plausible region in param- 
eter space has been rejected, although a large allowed 
range remains. Although we did not detect any signal 
in these 15 hours of observations, future detections may 
st rongly constrain the th eory. As found in Figs. 3 & 4 
in iBerezhko et all (|2006f) . the sensitivity of H.E.S.S. is 
much lower than the theoretical predictions. The future 
GLAST mission will also enable the lower energy range 
to be probed. Future large Cerenkov telescope arrays, 
such as CTA 19 will allow even more sensitive observa- 
tions to be made. 

7. CONCLUSION 

TeV gamma-ray observations toward the 400 year old 
remnant of Kepler's SN were made in 2005 April. Al- 
though a measurable flux of TeV-gamma rays had been 
predicted, we did not observe any statistically significant 
excess in this region, and the constraints on the allowed 
parameter range have been discussed. Although a re- 
gion of parameter range has been rejected, more sensi- 
tive measurements are required to constrain the models 
further. 
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